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Abstract 

We make an evaluation of the lepton flavor violating (LFV) decays r — > ^(1^*^, -P"), 
where £ = e or /i and V^{P^) is a neutral vector (pseudo-scalar) meson, in the context 
of unparticle physics. The constraints are investigated systematically on the related 
coupling parameters from all the available experimental data, and the parameter values 
are specified appropriately. The results show that whereas over the whole parameter 
space allowed by experiments all the r — > iP^ modes have a branching ratio too small 
to be measurable experimentally, in a large subspace as observed all of the r 
modes get simultaneously a branching ratio as high as 0(10"^'^ — 10~^), which is 
reachable at the LHC and super B factory. The important implications are drawn. 
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1 Introduction 



In the Standard Model(SM), massless neutrinos of different families are not mixed so that 
lepton flavors are made exactly conservative, or speaking, a lepton flavor violating (LFV) 
mode is forbidden absolutely in the SM. If the neutrino oscillation phenomenon takes place 
actually, we can affirm that the neutrinos are of a nonvanishing mass and thus lepton flavor 
conservation would be broken. Even so, LFV processes are still highly suppressed because 
of the smallness of neutrino masses. Hence, any distinct LFV signal can be deemed an 
indication of new physics beyond the SM. Recently, there has been an increasing interest in 
LFV physics. The current status of this subject is reviewed in [1]. 

Given the fact that the operators responsible for LFV transitions could be provided by 
most of the existing models beyond the SM, LFV phenomena could be explored in various 
theoretical frameworks. Most of efforts have been devoted to an investigation about LFV 
decays [SlEllllElEllTllHllQ] and lepton anomalous magnetic moments ((7 — 2) [H [TOl [11], [121 
[T3l [H]. A large branching ratio is predicted for r —>■ 3£, i'j and i{V^, P^) (where i = e or 
H and V^{P^) is a neutral vector (pseudo-scalar) meson) in some models such as the MSSM 
framework [3], SUSY seesaw mechanism [1], SUSY-GUT scenario [5] and type-Ill seesaw 
model [6]. Of all the existing discussions on LFV, those based on unparticle theory [15] are 
especially intriguing, because LFV processes can proceed at a tree level in this approach. 
The phenomenological implications of unparticle physics have been discussed intensely for 
the LFV transitions /i — > 3e [7], /i ^ 67 ^ and r ^ /i7 [9], electron and muon {g — 2) 
[SI [ini [m [121 [13] and collider physics [HI [16] . Besides, the effects of unparticle have been 
explored on hadronic processes [121 [13 [HI [IHl [20l [21] . More interestingly, the experimental 
constraints have been investigated on some of the unparticle coupling strengths and the 
important results have been obtained. For a mini-review on unparticle phenomenology one 
can be referred to [22] • The recent progress in unparticle physics can be found in [23]- On 
the other hand, a continuous experimental search has already performed for various LFV r 
decays. Very recently, an updated measurement has been reported on Brij — * IV^) [21] and 
Br{T — > a.K'^) [2S]- The estimated experimental upper limits on the branching ratios are in 
the range a few x(10~^ - 10"^) at 90% confidence level, for r £7 and ^{V^,P^) [26]. 

Though no clear signal has been detected in the current extensive search for LFV decays, it 
is expected that the future LHC will probe r ^ 3yU and liV^ , P°) down to the 10~^ level, 
while a sensitivity of 10"^'' — 10~^ will be reachable for a search for r 3^, £7 and £(K°, P°) 
at the super B factor [27] . 

Motivated by the recent progress in unparticle phenomenology and good prospect of 
the experiments on LFV r decays, in this Letter we intend to make an assessment of r 
^{y^^ P^) in the context of unparticle physics to understand the possibility to discover them 
in the future experimental searches. 
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This Letter is organized as follows. In the following section, on a brief introduction of the 
basic concepts of unparticle physics, we address the effective models we use for describing 
unparticle interactions with the SM particles and make a simple discussion. Section 3 is 
devoted to a derivation of decay rates for r i{V^,P^). A detailed parameter discussion 
and numerical evaluation is presented in section 4. The final section is reserved for summary. 



2 Effective Interactions 

The scale invariance in the conformal field theory, although not an exact symmetry of 
nature, might play an important role in exploring new physics beyond the SM. It prohibits 
strictly any particles with a definite nonzero mass from manifesting themselves and thus is 
broken in the SM. But there could be a sector, which is exactly scale invariant and interacts 
very weakly with SM particles at a scale much beyond the SM one. On the basis of a 
previous study [28], Georgi [15] suggests that there exist, in a very high energy theory, SM 
fields and BZ fields with a nontrivial infrared fixed point. These two sectors interact with 
each other by exchanging particles with a very large mass Mu. Below this mass scale, the 
heavy particles can be integrated out, resulting in the following local interactions: 

^ OsmObz, (1) 



j^dsM+dBZ—i 



where Osm is a SM operator with mass dimension dsM and Obz an operator with mass 
dimension dsz built out of BZ fields. When the energy scale runs down to a certain scale 
Au, at which the scale invariance in the BZ sector emerges, the renormalizable couplings 
of the BZ fields bring about dimensional transmutation. Then below this scale the BZ 
operators match onto unparticle ones and correspondingly, the interactions in (1) match 
onto an effective interaction of the form 



\dBZ-du 
^dsM+dBZ—'i 



OsmOu, (2) 



with Cu being a coupling coefficient and du the nonintegral number scale dimension of the 
unparticle operator On. 

Scale invariant unparticle stuff bears the characters strikingly other than those of ordi- 
nary particles. In particular, scale invariance can be used to fix the two-point functions of 
unparticle operators and further their propagators. The resulting propagators read, 

/A p^i pu 

d'xe^^-{0 I T^xWrn I 0) = ^^—^^(-9'" + -^)i-P' - ^e)'"-', (3) 

for a transverse vector unparticle, and 

d'xe'^'-^O I T[0u{x)0um \ 0) = ^—^^(-P' - lef--^ (4) 
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for a scalar unparticle. The coefficient Adj^ is given by 

167r5/2 T{du + l 



Ad,. 



(5) 



Since the matching procedure from the BZ operators to unparticle ones is unknown, 
unparticles may interact with SM particles in many possible ways. In the present case, we 
would like to use the effective coupling forms suggested by Georgi [15] . Then the interactions 
of a vector unparticle with the charged leptons can be expressed uniformly as 



u 



^ ^p.p. Ag/j- 



A, IP A 



A. 



Aru 



A. 



7 



(6) 



where a left-hand lepton vector is introduced, and all the related coupling constants 
Xij are arranged in a 3 x 3 matrix Ve and are treated as a real number. These coupling 
constants are in general viewed as a free parameter. A hierarchical relation, however we may 
conceive, does exist among some of them, because the LFV operators might be suppressed to 
a different degree by a small factor. We postulate that the following relations are respected: 
> Xr^i > Are aud A^^ > A^e- In fact, such relations could be accommodated by the 
existing experimental data, as will be seen later. 

The unparticle interactions with quarks could be discussed in parallel. Since only three 
light quarks are involved in the present situation, it suffices that we confine ourself to the 
former two generations. We have 



C 



u 



2kl[''^ULi^yuUL0i, 



2A 



u 



ul, cl 7^ 



Aun A, 
Aeii A, 




(7) 




2A^-'^"( dL, SL 



which describe the unparticle interactions with up-and down-type quarks, respectively. For 
the related flavor conserving couplings Agg (g = u, d, s) and flavor changing one \sdi we 
assume them to comply with the numerical relationship Au„ ~ A^^ ~ Xgs > A^^. 
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Correspondingly, the effective interactions involving scalar unparticle are of the following 
forms: 



2A 



u 



/ ^'pp -^e^i -^pr \ 



'lie 



XL A' , A' 



\ Ke Kfi Kt / 



\ 








/ 





(9) 



U 



2A^""hiL, CL 




(10) 



D 



7/. 



a: 



Al 




'u, 



:ir 



where it can likewise be assumed that there are the coupling hierarchies which are of the 
same structures as the corresponding ones suggested in the vector unparticle cases, and it 
should be understood that the scale dimensions have been set identical for the two different 
types of unparticles. 



3 Calculation of Decay Rates 

Now we embark upon calculating the decay rates for r ^{V^, P^) with the effective 
interactions (6) — (11). It is easily noticed that the scalar (vector) unparticle does not couple 
with a single vector (pseudo-scalar) meson. Then the decays r — > iV'^ proceed via just the 
vector unparticle, while the r —>■ iP^ transitions do by only the scalar unparticle. 

In order to discuss the vector unparticle mediated decays r — ^ IV^, we could take t ^ ficj) 
as an illustrative example. From the Feynman diagram plotted in Fig. [H we can write down 
the transition amplitude as 




Figure 1 : Feynman diagram of the vector unparticle- induced LFV decay r — > 
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A{t fJicj)) 



A. 



X 



2 sin((iev7r) 
|s7,(l- 7^)510), 



p2 



(12) 



where is the four momentum of the unparticle, and we have employed the ideal mixing 
scheme for the u — (p system. Further, the above expression can be simplified as, 



A{t — ^ 11(f)) 



-m, 



le 



(13) 



using the standard definition (0|s7^s|O) = m^f^e*^, with m^, and e being the correspond- 
ing mass, decay constant and polarization vector, respectively. After summing over the spins 
of the final states and averaging over the spins of the initial state, the decay width is derived 
as 

\P\ 



r r 



5^|^(r^/i0)P 



(14) 



spin 



where p stands for the momentum of the outgoing particles in the r rest frame, and 

n 2 



5^|^(r^/i0)P 



spin 



A2('iM-i) sin(dw7r) 



m^f^\{-ml-ie 



\du-1\ 



ml 



+ ml - 2m| 



(15) 



Using (13) — (15) and making a simple algebraic manipulation, the decay rates for the 
other r — > £V^ modes are easily achieved. Here we do not give them any more. 

For the scalar unparticle mediated decays r iP^, the hadronic matrix elements 
(P°(p) 1^17^75^2 1 0) enter into the expressions for the decay amplitudes. As usual, in the 
tt" and K^lK'^) case these matrix elements are parameterized, for instance, as 



{K'{p)\d^,j,s\0)=tfKP, 



(16) 



with Pfj^ and fx, respectively, being the four momentum and decay constant of the meson. 
In contrast, the t] and t]' situation is much more complicated because of mixing. The relevant 
decay constants are defined by 



(M(p)|g7'^75g|0) = {M{p)\sYl5s\0) = ift.p'', 



(17) 



where M = rj or rj' and q = u oi d. We would like to consider the rj — t]' mixing effect in the 
Feldmann-Kroll-Stech (FKS) scheme [29]. In this scheme the physical meson states \ri) and 



\ri'), in term of the parton Fock states |?7g) 
as 



1^) 



COSi 

sin ( 



\uu + dd)/\/2 and |?7s) = \ss), are decomposed 

-sin0 \ / \r]g) 
COS0 j \ \r]s) 
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where (j) is the mixing angle. Furthermore, by defining the two basic decay constants fq and 
fs as 



{Vq{p)\Ql^l5q\0) = -^fqP^, {Vs{p)\si^l5s\0) = ifsP^, 



(19) 



we have the following relations: 



frj = fq COS 0, = - fs sin 

= fq sin (f), f^f = fs COS 0. 



(20) 



With the aid of the data fitting results fjf^ = 1.07, fs/fn = 1-34 and = 39.3° ± 1.0°, the 
desired values of the decay constants /^,, and f^, can be achieved f2U\ . 

At present, the decay rates for r IP^ could be calculated with the known decay 
constants. As in the vector meson case, we illustrate our findings of r ^ IP^ with the 
resulting expression for the decay width in the r ^K^ case. 



r(r ^ /iJ-sT" 



\P\ 



0m2 



spin 



with 



^T/i^sd -^du n 2 I 2 

-fRrriKi-mj^ 



le 



2h?du sin((iz^7r) 
X [m^Jiil + 7^)r + m^7l(l - 7^)r] 



and 



0m2 



spm 



n 2 



A2rfM sin(rfz^7r 



-fKml\{-ml-teY--^\ 



X 



/ 2 I 2\f 2 I 2 



2 2 



(21) 



(22) 



(23) 



4 Parameter Discussions and Numerical Evaluations 

For a quantitative analysis of these LFV decays, we need to make a detailed discussion 
about the various parameters involved in the calculation. 

The mass parameters associated with the present calculation have been well known, the 
decay constants of the related light mesons have been determined better experimentally too. 
All those are listed in TabiU 

The non-integral scale dimension du is calculable in principles, but difficult to estimate 
in practice. However, it might be limited to 1 < rf^^ < 2, which is to be used here, from 
the unitarity [IS, [30] and convergence condition. As concerns the scale parameter An, we 
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could let it range from 1 TeV to a few TeV, because ones expect generally that a certain 
new physics, if it exists, should appear at such energy region. 

Our main concern, of course, is how the underlying coupling constants take their values. 
For having a knowledge about the couplings A^^ and A^e? ones have to make a correlation 
discussion for the various LFV processes where these couplings are involved. Unfortunately, 
the currently available experimental data are not sufficient to provide them with a decisive 
parameter space. For the relevant unparticle-quark couplings, though we can extract them 
from the experimental measurements on some hadronic processes in a certain data fitting 
way, the uncertainties, among other things, in the long distant QCD parameters would affect 
greatly the accuracy of extraction. 

As the case stands, it is needed to work at a level of order of magnitude, as we make 
a choice of parameter sets from the regions allowed experimentally. Before starting our 
discussion, a couple of explanations are in order: (1) Both vector and scalar unparticles 
could in general be responsible for a LFV transition. Including simultaneously contributions 
of both the unparticles can make not only the results have a large uncertainty but also 
the calculation extremely complicated. In the following parameter discussion we assume 
that they two contributes separately, as done in many studies, and consider only the vector 
unparticle cases. Also, we suppose that the corresponding coupling strengths are the same for 
the scalar and vector unparticle interactions. (2) We know that the unparticle parameters of 
an effective interaction contains Kui du and the coupling constant A. The resulting transition 
amplitude for a process depends on the parameter function fdi^{X,Ai(). With the function 
values extracted from an experiment, which are generally relevant to du, the values of A 
can be determined at any du, the results being, of course, dependent on Ak. Accordingly, 
these coupling values, though changed with Ak, correspond to one and the same /(A,A^) 
of a fixed value. If we want to make a theoretical prediction with these extracted coupling 
parameters, we could work at an arbitrarily chosen Ak. The final results must have nothing 
to do with Ak, for the same /(A, A^^) enters, which keeps its value unchanged for different 
Ak- For convenience, we will work at Ak = 1 TeV. (3) The unparticle couplings Age and A^e, 
as two important inputs in our parameter determination, have been investigated in detail 
in [ini El] and [Hj, respectively. In the region 1.5 < dK < 2, the resulting bounds on Age 
and A^e are available for the present case. From the findings obtained by a study on the 
inviable positronium decays [31], we deduce easily that Age < 10^"^ for (Ik = 1-5, Age < 10~^ 
for dK = 1-6 and Age < 10^^ for dK > 1-7. Moreover, A^g has a negligibly small number, 
as required by the experiments on the /i — e conversion in heavy nuclei [8] , so that we can 
set it to zero. These constraint conditions will be used below to restrict other unparticle 
couplings. In the region 1 < dK < 1-5, the study indicates that a more stringent restriction 
on Aee comes from the precise measurement on long-ranged spin-spin interaction of electrons 
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[3T] . However, the results are not directly applicable and a revaluation is needed. In reality, if 
we work in the present context we have to assess not only Age but also the related unparticle- 
quark couplings in the region 1 < < 1.5 in which these unparticle parameters are less 
known. It is possible to make such an assessment, however goes beyond the scope of this 
work. We will choose 1.5 < rf^^ < 2 as our work region. 

We investigate, to begin with, the possible regions of Xr^ and A^e allowed by the existing 
experiments [26]. The authors of [9] manage to understand the parameter region of A^^ in 
a scalar unparticle model by a combined analysis of the muon g — 2, r fi'y and r 3fi. 
They find that as du > 1.6 the muon g — 2 experiment demands that A^e have a negligibly 
small number, which is in agreement with what is required by the experiments of the fi — e 
conversion in heavy nuclei [8], and at least one of A^,-,^(/i be of O{10~^ — 1). Including 
further the possible constraints from the experimental data i?r(r 3fi) < 3.2 x 10~^ and 
Br{T fi'-f) < 6.8 X 10~^, they conclude that one of the two couplings is of 0(10^^ — 1), 
while the other is at or below order 10"^. These constraints are possibly weak, because 
in the derivation the fi loop is assumed to dominate in the r — > /i7 transition so that the 
contribution of the virtual r particle is not included. However, the same conclusion can yet 
be drawn in disregard of the constraint of r ^ fi'-f. We make the same investigation within 
the present framework by means of the experimental observations of the muon — 2, r — > 3/i 
and r — > fie^e~ (with Brij fie^e~) < 2.7 x 10~^), and find that A,-^ can range from 
O{10~^) to O{10~^) if A^^ takes a larger value of O{10~^ — 1), and vice versa. From these 
possible parameter regions we can pick out our preferred parameter sets: (1) A,-^ = 10~^ 
and A^^ = 10"^, for du = 1.5. (2) A^^ = IQ-^ and A^^ = 10~\ for du = 1.6. (3) A^^ = 10"^ 
and A^^ = 1, for du > 1.6. At this point, we must emphasize the fact that the current 
experimental data on the tau g — 2 [26] do not provide more about the couplings involving r 
lepton than we get above and below, because of the existing sizable uncertainty which allows 
us to do theoretical calculation within a considerably large space of parameter. 

As far as A^e is concerned, the parameter regions allowed by /i ^ 67 have been evaluated 
in [8], However, a consistent evaluation requires us to consider a combined constraint from 
the processes /i —>■ 67, r ^ 3e and r efi~^fi~ as well as electron g — 2. From the 
experimental measurements Br{fi 67) < 1.2 x 10^^^, i?r(r 3e) < 3.6 x 10~® and 
Br{T efi^fi^) < 3.7 x 10^^, it follows that A^e can be limited to the region A,-e < (9(10~^) 
for du = 1.5 and 1.6, while the resulting upper limits can basically remain at order of 
10~^ for du > 1.6. The constraint is achievable from the electron g — 2 experiment too, 
by making a replacement of the corresponding parameters in the expression for the muon 
g — 2 and then confronting the result with the numerical deviation between the SM estimate 
and experimental measurement |Aa;| < 15 x 10^^^ [32]. But no new results are found. In 
the numerical evaluation we will use A^e = 10~^ for du = 1.5 and 1.6, and A^e = 10"'^ for 
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du > 1.6. 

In passing, it is attractive to examine the possible region for A^r using the experimental 
bounds Br{T fi'-f) < 6.8 x 10~* and i?r(r 67) < 1.1 x 10"'', along with the various 
constraint conditions obtained already. The results show that t ^ fi'-f furnishes a stronger 
restriction A,-,- < (9(10) as du > 1.5 and therefore the possibility of a sizable unparticle-tau 
coupling strength cannot be ruled out. Then we can conclude, according to the present 
study, that our hierarchy assumptions A,-^ > A,-^ > \re and A^^ > A^e are acceptable at 
least for the existing LFV experiments. It remains to be seen whether such relationships 
are true or not. We can believe that the future precision measurement on the tau g — 2 [33] 
would help to clarify this issue. 

To turn to the discussion about the unparticle-quark couplings. The existing constraints 
on them come mainly from the studies on some inclusive [2D] and exclusive [T5| [TU] decays 
of B mesons and neutral meson mixing systems [121 [13, [ISl EH]. One expects that the 
inclusive process B would provide a stringent constraint on new physics effects, as 

a result of the good agreement between the experimental measurement and SM prediction 
on Br{B ^sl)- However, it is not always this case in the face of unparticle effects 
|20] . The sensitivity of Br{B Xg'y) to the coupling parameters weakens as du > 1.5. 
As it is, the constraints would become considerably weak as du > 1.7 so that a sizable 
unparticle-quark coupling strength is allowed. To have more understanding of the unparticle 
parameters, in [19] the impacts of unparticle are analyzed on B^^s — Bd^s mixing processes 
and exclusive channels B tttt, vri^ , and especially a detailed data fitting is carried 
out for the B ^ nir, irK with the constraints of B^^s — Bd^s mixing. The fitting results 
with du = 1.5 demonstrate that there is a large coupling strength of 0(10"") (n = 0, 1) for 
the flavor conserving interactions, which is compatible with the findings in the B — > Xg'^ 
situation. What is particularly interesting is that A„„ and \dd turn out to be at the same 
order of magnitude and a relatively small number is implied for Ag^, as expected by us. 
It is claimed that with the yielded optimized values of parameters, the existing various 
discrepancies may be explained between the SM predictions and experimental data. We 
think that these constraints, though subject to an estimate of uncertainty, could serve as a 
valuable reference for us to select proper parameter values. We assign the following numbers 
to the related couplings: A„„ = A^^ = A^s = A ~ 10"" {n = —1,0,1,2) and A^^ = 10~"A 
{n = 1,2), for du = 1.5. In the region du > 1.5, little is known about them. Nevertheless, 
numerous studies show that when du increases, the ranges allowed experimentally become 
large for unparticle-lepton couplings. The same should be true of the quark case, for the 
coupling forms are the same in the two situations. Taking this point into consideration 
and for simplicity, we suggest that these unparticle-quark couplings remain unchanged in 
the region du > 1.5, due to a certain stringent restriction condition. This is equivalent to 
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Table 1: Summary of the leptonic and hadronic parameters (in units of MeV). 





fK 


fp 






f<f> 


/. 
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130 


160 


209 


195 


217 


231 


139 


174 




rriK 


nip 




TTlK* 








130 


498 


770 


782 


892 


1020 


547 


958 








TJlr 
















nil 


105 









a conservative estimate. In addition, in the case of r — £(p'^,7r°) we need to confront a 
combination of two couplings \uu — ^dd- Since we are discussing the unparticle couplings at 
a level of order of magnitude, it is sound to set A^^ — A^^ at the same order as Xuu,dd- 

Now we are a position to make a numerical evaluation. In the first place, we can notice 
that for both r iV^ and r iP^ an approximate order of magnitude relation exists 
between the branching ratios, with our selected coupling parameters. In the r iV^ 
situation, we have the following observation: 

Br{r — s> fip^) ~ Br{r jiuj) ~ Brij ficj)) 
> Br{T ep^) ~ i?r(r euj) ~ Br{T — > ecf)) 
> Br{r pK*°{K*^)) > Br{T eK*\K*°)), (24) 

if neglecting the mass difference between muon and electron and SU (3) breaking effects in 
the hadron parameters. A similar relation holds approximately for r iP^. However, it 
would suffer from a large SU{3) breaking correction. The numerical calculations denote 
that these order of magnitude relations are, indeed, respected better for r £V^ than for 

Let us take a closer look at the behaviors of the branching ratios in the parameter spaces 
we adopt. It is clearly seen that the parameter region A < 1 is allowed by the experiments, 
while the region A > 1, where the branching ratios for all the r fiV^ go beyond their 
experimental upper limits, is prohibited. The allowable parameter sets are fixed as: (I) 
A = 10~2, 1.5 < du <2; (II) A = 10~\ 1.55 < du < 2; (III) A = 1, 1.85 < du < 2. Over 
these parameter areas all the r — > iP^ modes show a branching ratio less than (9(10~^°), 
which is far from the experimental reach. We will focus our discussion on the r iV^ 
case. For the set I, the branching ratios are of orders 10~^^ — 10~^ for r pV^, compared 
with the numerical region for Br{T eV^) 10~^^ — 10~^^. In the set II case, whereas the 
r fiV^ modes have a branching ratio ranging from 10"^^ to 10~^, the numerical results 
for Br{T eV^) are located between 10~^^ and 10~^. If the set III is used, the numerical 
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values for i?r(r f^V^) vary from lO"^*^ to 10~^, while those for i?r(r eV^) do between 
10~^^ and 10~^°. To illustrate the dependence of i?r(r iV^) on A and du, we can typically 
consider the r — ^ fi(j) case in which the behaviors are shown of i?r(r fi(f)) in some parameter 
regions in Fig. 2. Albeit the branching ratios turn out to be sensitive to du and A, there is 
still a large parameter region, as will be seen, in which for any du almost all the r — > fiV^ 
modes have a branching ratio as large as (9(10~^° — 10"^), which are expected to be reachable 
at the LHC and super B factor. 

The parameter regions experimentally favorite can be summarized as: (I) A = 10~^, 
1.5 < du < 1.8; (II) A = 10~\ 1.55 < du < 1.9; (III) A = 1, 1.85 < du < 2. In these 
regions of the allowable parameter space, as a matter of fact, all the r fiV^ modes except 
r — > fiK*^{K*^) are accessible experimentally, and in some subregions the same observations 
can be obtained for r fj,K*'^{K*^) and r — * e(p°,u;,0). Only r eK*^{K*^) exhibits a 
branching ratio below (9(10~^°). The partial findings from these parameter regions, together 
with the current experimental upper limits on them, are collected in Tab. 2. 

Table 2: Some selected numerical results for i?r(r — > ^V^). The corresponding pa- 
rameter sets (A^^, Ke, A, \sd) are (lO'^, 10"^, 10"^ lO'^), (IQ-^, IQ-^, lO-\ lO'^), 
(10-2, 10-^ 10"\ 10-2) and (10~2, 10-^ 1, 10"^), respectively, for du = 1.5, 1.6, 1.7 and 
1.8, and 1.9. 



du 

Mode\^ 


1.5 


1.6 


1.7 


1.8 


1.9 


EXP.UL 


T fip^ 

T jlUJ 
T — S> /i</) 

r ^ /iA'*° 
r ^ nK*^ 


8.0 X 10-9 
6.9 X 10-9 

8.1 X 10-9 

1.2 X 10-^° 
1.2 X 10"^° 


2.6 X 10-s 
2.2 X 10-^ 
2.9 X 10-* 
4.2 X 10-1° 
4.2 X 10-1° 


6.3 X 10-« 

8.1 X 10-« 

1.2 X 10-^ 
1.6 X 10-9 
1.6 X 10-9 


4.3 X 10-9 

3.8 X 10-9 
6.1 X 10-9 

7.9 X 10-11 
7.9 X 10-11 


3.6 X 10-s 
3.2 X 10-* 

5.7 X 10-* 
6.9 X 10-1° 
6.9 X 10-1° 


6.8 X 10"* 

8.9 X 10"* 
1.3 X 10-^ 
5.9 X 10"* 
1.0 X 10"^ 


T ep^ 
r — i> ecu 
r ^ e0 
T eK*^ 
T ei?*° 


8.1 X 10"" 
7.0 X 10-" 
8.3 X 10-" 

1.2 X 10-12 
1.2 X 10-12 


2.6 X 10-1° 
2.2 X 10-1° 
3.0 X 10-1° 
4.2 X 10-12 
4.2 X 10-12 


9.4 X 10-1° 
8.2 X 10-1° 
1.2 X 10-9 
1.6 X 10-11 
1.6 X 10-11 


4.4 X 10-11 

3.8 X 10-11 
6.2 X 10-11 

7.9 X 10-13 
7.9 X 10-13 


3.6 X 10-1° 
3.2 X 10-1° 
5.8 X 10-1° 
7.0 X 10-12 
7.0 X 10-12 


6.3 X 10"* 
1.1 X 10-^ 
7.3 X 10-* 
7.8 X 10"* 
7.7 X 10-* 



So far all the numerical calculations are performed with the fixed coupling values. How- 
ever, these parameter values, as has been emphasized, should are understood as an order of 
magnitude and thus we have to consider the effects resulting from the variations of coupling 
parameters within their individual orders of magnitude. We have examined such effects. 
With the fixed values of A^^ and A^e, A dependence of Br{T — > /i0) is plotted in Fig.3. 
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1.80 
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Figure 2: du dependence of Br(T — > ^(j)) with the different A values. The horizontal line denotes 
the present experimental upper bound on i?r(r — > ^(j)). 




Figure 3: A dependence of i?r(r — > fi4>) with the different dy. The horizontal line denotes the 
present experimental upper bound on i?r(r 
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Obviously, the numerical results can change by up to two orders of magnitude, when du 
remains fixed and A ranges from 0.01 to 0.1. The similar situation appears for 0.1 < A < 1. 
Exploring further the case where all the related parameters vary simultaneously, we would 
have a much larger numerical range. But this also indicates that there could are more theo- 
retical results which are within the experimental reach. Those listed in Tab. 2 are only some 
estimated lower bounds on i?r(r — ^ IV^). In fact, when the couplings A,-^ and A^e change at 
the same time within their respective ranges, we can work within the expanded parameter 
regions for A and du- (I) 0.01 < A < 0.04, du = 1.5; (II) 0.01 < A < 0.1, 1.5 < du < 2; 
(III) A = 0.1, 1.55 <du<2; (IV) 0.1 < A < 1, 1.75 < du < 2; (V) A = 1, 1.85 < du < 2. 
In these parameter spaces there are many favorable subspaces as observed, in which all the 
r — > iV^ modes can get simultaneously a branching ratio of (9(10^^° — 10~*). By contrast, 
all the r iP^ modes remain still inaccessible to experiments over these areas. 

Conversely, the restrictions can be inspected on the coupling parameters from the exper- 
iment data on 5r(r ^(^°, P^))- By limiting ourself to the parameter area 0.01 < A < 1, 
the constraints on At-^ and X^e are observed to be weaker in a large region of the parameter 
spaces for A and du, and rather loose in some subspaces like A = 0.01 and 1.6 < du < 2, 
when compared with those presented above. 

Again we stress that all the presented numerical results, despite achieved at A^^ = 1 TeV, 
maintain unchanged as Au increases, and besides, the scalar unparticle mediated r — > iP^ are 
evaluated using the same scale dimensions and coupling strengths as in the vector unparticle 
case. Even if we regard these corresponding parameters as independent of each other, it 
is yet difficult to get a interesting result for r — iP^, since it is hardly conceivable that 
the related coupling constants have a considerably sizable number in such a case. Once 
the coupling parameters become better understood in the whole du region 1 < rf^^ < 2, we 
could make a more complete and reliable assessment of these LEV processes. But it seems 
likely that with the parameter sets specified adequately the hierarchical relation i?r(r 
£V^) ^ Br{T iP^) will be kept valid, although the branching ratios alter with change in 
parameter values. Of course, to do calculation with different unparticle coupling scenarios 
would in general lead to different results. It is desirable to enquire into these LEV processes 
in other unparticle coupling schemes. 

Our findings for r — > iV^ appear to be comparable with some of the existing estimates 
[31 m [5], [6] . Nevertheless, in the r iP^ case we have a branching ratio much less than those 
for T — ^ iV^, presenting a striking contrast to the predications of the other approaches. 
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5 Summary 



We have made a detailed analysis for the unparticle induced LFV decays r — > i(y^, P°) 
in an effective model with a hierarchical relation suggested among some of the coupling 
constants. 

To get a consistent and believable assessment, all the available experimental data have 
been used to constrain the unparticle couplings. From the obtained constraint conditions, 
the parameter values for the related couplings have been specified appropriately. As a by- 
product, it is found that a sizable Xrr is allowed by the current experimental data, and 
our hierarchy hypotheses A^r > > Ke and A^^ > A^e can be accommodated by these 
constraint conditions. 

We have evaluated the branching ratios and examined the possibility to experimentally 
discover these modes in the near future. In the parameter region A > 1, all the r — > fiV^ 
modes have a branching ratio exceeding their individual experimental upper limits. The 
experimentally allowed regions for A and are determined approximately as: (I) 0.01 < 
A < 0.04, du = 1.5; (II) 0.01 < A < 0.1, 1.5 < du < 2; (III) A = 0.1, 1.55 < du < 2; 
(IV) 0.1 < A < 1, 1.75 < du < 2; (V) A = 1, 1.85 < du < 2. In many regions of these 
parameter spaces, for all the r — > iV^ modes we can have simultaneously a branching ratio 
of orders 10"^'^ — 10"^, which are expected to be accessible at the LHC and super B factory. 
Compared with the r — iV^ case, all the r — > modes show a branching ratio beyond 
the experimental reach. 

Also, we have inspected the limits imposed on the couplings A^/^ and A^e by the experi- 
ments on r ^ i{V^,P^), observing that there is a looser bound than those yielded by the 
other available LFV experiments, in a large subspace of the parameter spaces 0.01 < A < 1 
and 1.5 < du < 2. 

It is explicitly too early to draw a final conclusion whether these LFV decays are ob- 
servable experimentally. We have to await the improvement in experiment and progress in 
unparticle phenomenology. Different from the predictions of the other new physics models, 
however, the unparticle approach gives the numerical relation i?r(r — iV^) ^ i?r(r — > iP^), 
with the implication that there is a greater discovery potential of r — > £V^ than that of 
r —>■ iP^ in future experiments. If this gets confirmed in the future experimental searches, 
the present research is perhaps instructive in identifying whether or not these LFV processes 
are induced or dominated by unparticles. 
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